Where is the jet quenching in Pb + Pb collisions at 158 AGeV? 
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Because of the rapidly falling particle spectrum at large px from jet fragmentation at the CERN 
SPS energy, the high-pT hadron distribution should be highly sensitive to parton energy loss inside 
a dense medium as predicted by recent perturbative QCD (pQCD) studies. A careful analysis of 
recent data from CERN SPS experiments via pQCD calculation shows little evidence of energy 
loss. This implies that either the life-time of the dense partonic matter is very short or one has 
to re-think about the problem of parton energy loss in dense matter. The hadronic matter does 
not seem to cause jet quenching in Pb + Pb collisions at the CERN SPS. High-pr two particle 
correlation in the azimuthal angle is proposed to further clarify this issue. 



Hard processes have been considered good probes 
of the dense matter which is produced in high-energy 
heavy-ion collisions and is expected to be in the form of 
deconfined quarks and gluons or a quark-gluon plasma 
(QGP) at high energy densities. These processes hap- 
pen in the earliest stage of the collisions and therefore 
can probe the properties of the dense matter in its early 
form, whether a QGP or not. Furthermore, their pro- 
duction rates can be calculated with reasonable accuracy 
within pQCD parton model and has been tested exten- 
sively against vast experimental data in p + p and p + A 
collisions. These calculations Q incorporating minimum 
amount of normal nuclear effects (nuclear modification 
of parton distributions Q and Cronin effect ||) then 
provide a clean and reliable baseline against which one 
can extract signals of the dense matter. In this paper, 
we investigate what high-py particles from jet fragmen- 
tation tell us about the dense matter formed in Pb + Pb 
collisions at the CERN SPS. 

Like other hard processes, large transverse momen- 
tum parton jets are produced in the early stage of high- 
energy heavy-ion collisions. They often have to travel 
through the dense matter produced in the collisions and 
finally hadronize into high-p^ particles in the central 
rapidity region. Recent theoretical studies show 
that a fast parton will lose a significant amount of en- 
ergy via induced pQCD radiation when it propagates 
through a dense partonic matter where the so-called 
Landau-Pomeranchuk-Migdal coherence effect becomes 
important. If this picture of parton energy loss can be 
applied to large transverse momentum parton jets in the 
central rapidity region of high-energy central A +A colli- 
sions, one should expect a leading parton to lose energy 
when it propagates through a long-lived dense matter. 
Since the radiated gluons will eventually become inco- 
herent from the leading parton which will fragment into 



large-pT hadrons, one then should expect a reduction of 
the leading hadron's pt or a suppression of the large- 
p T particle spectrum J§-[h}. At the CERN SPS energy, 
high-pT jet or particle production (jpr > 3 GeV/c) is 
very rare and the power-law-like spectrum is very steep 
because of the limited phase space. It should be espe- 
cially sensitive to any finite energy loss. 

The single inclusive particle spectrum at large pt in 
high-energy p+p or p+p collisions can be calculated in a 
pQCD parton model with the information of parton dis- 
tributions jlj] and jet fragmentation functions Jj5fl from 
deep-inelastic e + p and e + e~ experiments. This is one 
of the early successes of the QCD parton model |ll] [l3| ]. 
It was already pointed out that the initial transverse 
momentum before the hard scattering is very important 
to take into account at lower energies and can signifi- 
cantly increase the single inclusive differential cross sec- 
tion. The initial parton transverse momentum can be 
studied in detail via Drcll-Yan (DY) |@-|§, 7+jet and 
7 + 7 production in p + p collisions. 

To the lowest order of pQCD, the single inclusive par- 
ticle production cross section can be written as [13|, 



dycPpr 



J dx a dx h d 2 k aT d 2 k bT g p (k aT ,Q 2 ) 

abed 



g P (k b T, Q 2 )fa/ P (x a , Q 2 )fb/ P (xb, Q 2 ) 

D° h/c (z c ,Q 2 ) da , 



dt 



(ab — > cd), 



(1) 



where x a $ are the fractional energies and fc a ,t>T the 
initial transverse momenta of the colliding partons. 
da/dt(ab —> cd) are the differential elementary parton- 
parton cross sections |D|]. K w 2 is used to account 
for higher order corrections p(J and Q = P c t — Pt/z c . 
We will use MRSD— ' parameterization for the par- 
ton distributions f a /p(x,Q 2 ) and BKK parameteriza- 
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tion for the jet fragmentation functions D^^ c (z, Q 2 )- We 
will use a Gaussian form for the initial-fc-r distribution 



l/(7r(fcj.) p ) exp(— fcj./(fcj,} p ) with a vari- 



ance (fc^)p = 1 (GeV 2 /c 2 ) + 0.2Q 2 a s (Q 2 ), where the 
Q-dependence accounts for initial kx from initial-state 
radiation ( or higher order 2 — > 2 + n processes) p|. 
The parameters are chosen to best fit the experimental 
data of high-p-r particle spectra at all energies Q . Be- 
cause of the introduction of initial parton kx, one of the 
Mandelstam variables for the elementary parton-parton 
scattering processes could vanish and cause the differ- 
ential parton cross sections to diverge in certain phase 
space points. We use an effective parton mass fi = 0.8 
GeV to regulate the divergence as in the early studies 
[p~T]| . The resultant spectrum is sensitive to the value of 
/i only at around pt ~ (J., where pQCD calculation is not 
reliable in any case. 

Shown in Fig. [j] is an example of the calculated 7r ± 
spectra in p+p collisions at E\ a h = 200 GeV. The agree- 
ment with experimental data is very good not only for 
the overall inclusive cross section but also for the iso- 
spin dependence as shown by the p-r-dependence of the 
7r~ /n + ratio in the inserted figure. Similar analyses have 
been carried out at other energies up to Fermilab Teva- 
tron The initial kx is less important and becomes 
almost negligible for the single inclusive parton spectra 
at these collider energies. 




1 



7 



3 4 5 6 
p T (GeV/c) 

Single- inclusive pion spectra in p + p collisions 
200 GeV. The solid lines are from pQCD calcu- 
lations and data from Ref. The inserted figure are the 
corresponding tv~ /tt + ratios. 
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In p + A collisions, there are two known nuclear ef- 
fects: nuclear modification of the parton distributions 
(EMC effect) [|| and nuclear enhancement of the large- 
Pt hadron spectra (Cronin effect) H. Both are caused 
by multiple initial scattering. We assume that the par- 



ton distributions per nucleon inside a nucleus at impact 
parameter b, 



fa/A(x,Q 2 ,b) = S a/A (x,b) 



fa/ P (x,Q 2 ) 



+ (l~j)fa/n(x,Q 2 ) 



(2) 



is factorizable into the parton distributions inside a 
normal nucleon and the nuclear modification factor, 
S a / A (x, b), for which we use the HIJING parameteriza- 
tion ||). This should be adequate at the CERN SPS 
energy where the dominant process at large pr is quark- 
quark scattering. 

One can explain the Cronin effect within a multiple 
parton scattering model J23|,^4]], in which the cancella- 
tion by the absorptive processes forces the nuclear en- 
hancement to disappear at large pr like l/p T and in the 
meantime causes a slight suppression of hadron spectra 
at small pt so that the integrated spectra do not change 
much . This allows us to take into account the effect of 
multiple scattering via a broadening of the initial trans- 
verse momentum, 



(k 2 T ) A {b) 



(4) P + K6)-1]A 2 



(3) 



where v(b) = a pp t A {b) is the average number of scat- 
tering the parton's parent nucleon has suffered and 
t A (b) is the nuclear thickness function normalized to 
J d 2 btA(b) = A. Since the Gaussian distribution is not 
a good approximation for the fcx-kick during the ini- 
tial multiple scattering, we found that we have to use 
a scale-dependent value, A 2 = 0.225 ln 2 (Q/GeV)/(I + 
ln(Q/GeV)) GeV 2 /c 2 , to best describe the available data 
from p + A collisions [^lj which allow about 10 — 20% 
uncertainty in the calculated spectra. For Q = 2 ~ 3 
GeV, A 2 = 0.064 - 0.129 GeV 2 /c 2 , which is consistent 
with the analyses of pr broadening for J/ty production 
inp + A 

Taking into account these nuclear effects which al- 
ready exist 'mp+A collisions, the single inclusive particle 
spectra in A + A collisions can be estimated as 



dyd 2 p T 



abed 

d 2 k aT d 2 k bT gA(kaT, Q 2 , r)g A (k bT , Q 
,Q 2 ,r)f b/A (x b ,Q 2 ,\b-r\) 
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d 2 b / d 2 rt A (r)t A {\h - r|) / dx a dx b 



|b-r|) 



fa/A(x a 

D° h/C (z c 



Q 2 )da 



■KZn 



At 



-{ah — > cd). 



(4) 



The initial-fcr distribution g A (kx,Q 2 ,b) is similar to 
that of a proton in Eq. (Q) with a broadened width given 
by Eq. (|^) which now depends on the impact-parameter 
b. 

For central A + A collisions, we limit the integration 
over the impact parameter to & max - Using the geometri- 
cal cross section of a hard-sphere nucleus, we determine 
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frmax by matching b 2 max /AirR 2 A (R A w 1.12A 1 / 3 fm) to the 
fractional cross section of the triggered central events in 
experiments. In Eq. ||, we actually use the Wood-Saxon 
distribution to calculate the thickness function tA(b). 

Shown in Fig. || are the calculated single-inclusive 
spectra for 7T° in central S + S (E lah = 200 GeV) and 
Pb + Pb (E lah = 158 GeV) collisions with (solid) and 
without (dashed) nuclear A^-broadening as compared to 
WA80 |3 and WA98 |§ data. Besides small effects of 
the nuclear modification of the parton distributions on 
the spectra at these energies, the dashed lines are simply 
the spectra in p + p collisions multiplied by the nuclear 
geometrical factor. It is clear that one has to include the 
fc-r-broadening due to the initial multiple scattering in 
order to describe the data. This is also consistent with 
the analysis by WA80 §7j. 
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FIG. 2. Single-inclusive n spectra in central S + S at 
Blab = 200 GeV and Pb + Pb collisions at £i ab = 158 GeV. 
The solid lines are pQCD calculations with initial-fcT broad- 
ening and dashed lines are without. The S + S data are 
from WA80 |3 and Pb + Pb data are from WA98 Q. The 
dot-dashed line is obtained from the solid line for Pb + Pb by 
shifting p T by 0.2 GeV/c. 

One can conclude from this analysis that the factor- 
ized pQCD parton model seems to work well for large-pT 
hadron production in A + A collisions. But one can also 
immediately realize that there is no evidence of parton 
energy loss as predicted by previous theoretical studies 
. If there is parton energy loss and the radiated glu- 
ons become incoherent from the leading parton, the ef- 
fective fragmentation functions should be modified such 
that the leading high-p^ particles should be suppressed 
as compared to p + p and p + A collisions J8|— p~0[ . At 
y = 0, parton energy loss can be directly translated into 
Pt reduction for the leading hadrons. To estimate the 
experimental constraints on parton energy loss, one can 



simply shift the pr values of the solid line for Pb + Pb in 
Fig. | by 0.2 GeV/c (dot-dashed line). Assuming 20% 
uncertainty of the calculated spectrum, one can quickly 
exclude a total energy loss AE < 0.1 GeV. With the 
transverse size of a Pb nucleus, this corresponds to an 
energy loss dE/dx < 0.02 GeV/fm. Detailed model cal- 
culations will give a more stringent limit pl| . This is in 
direct contradiction with the current theoretical studies 
of parton energy loss in dense matter and calls into ques- 
tion current models of energy loss. It also implies that 
there is not a dense partonic matter which exists long 
enough to cause parton energy loss. 

Most of the recent theoretical studies of energy 
loss are based on pQCD calculation for a single fast par- 
ton propagating through a large dense medium. If we 
assume that it is valid for a parton propagating through 
a deconfined medium, the absence of parton energy loss 
in the experimental data on high-pT particle spectra im- 
plies that either there is no such deconfined partonic 
matter being formed or it only lived for a very short 
period of time. Using the measured dEx/drj rs 405 GeV 
[ p9| in the central rapidity region of most central Pb+Pb 
collisions (%2 of the total inelastic cross section) one can 
estimate the initial energy density at tq = 1 fm/c to be 
about e = dE T / dr] / {-kt R\) « 2.9 GeV/fm 3 . This is 
an optimistic estimate assuming that the formation time 
of the dense matter is about 1 fm/c. Because of longi- 
tudinal expansion, the energy density will decrease like 
e/e = (r /r) Q . The value of a could range from 1 for 
free-streaming to 4/3 for hydro-expansion of an ideal gas 
of massless particles. Assuming a critical energy density 
of e c » 1 GeV/fm 3 , the system can only live above this 
critical density for about 2.2 ~ 2.9 fm/c. Equilibrating 
processes and transverse expansion certainly will reduce 
this life-time even further. During such a short time, 
a highly virtual parton has small interaction cross sec- 
tion before its virtually decreases through pQCD evo- 
lution. Therefore, a produced large px parton will not 
have much time to lose its energy before the dense mat- 
ter drops below the critical density. The recent theoret- 
ical studies are not applicable to such a short-lived 
system. Nevertheless, this analysis at least tells us that 
the life-time of the dense partonic matter must be short 
if it is ever formed in Pb + Pb collisions at 158 AGeV. 
Otherwise, it is difficult to reconcile the absence of par- 
ton energy loss with the strong parton interaction which 
drives the equilibration and maintains a long life-time of 
the initial parton system. 

One definite conclusion one can draw from this anal- 
ysis is that the hadronic matter in the later stage of 
heavy-ion collisions does not seem to cause parton en- 
ergy loss or jet quenching at the CERN SPS. This will 
make jet quenching an even better probe of long-lived 
initial partonic matter since it will not be affected by 
the hadronic phase of the matter. Because of its long 
formation time (r/ ~ 20 fm/c for a pion with px ~ 3 
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GeV/c), a high-py pion is only formed either after freeze- 
out or in a very dilute hadronic matter. Otherwise, in- 
elastic scattering with other soft pions can also cause the 
suppression of high-p^ particle spectra or apparent jet 
quenching. What is traveling through the hadronic mat- 
ter is thus a fragmenting parton whose interaction with 
a hadronic matter might be non-perturbative in nature. 
The pQCD estimate of parton energy loss is then not 
applicable here even though it might be adequate for a 
parton propagating in a hot QGP. The fact that a frag- 
menting parton does not lose much energy in hadronic 
matter might be related to the absence of parton energy 
loss to the quarks and anti-quarks prior to DY hard pro- 
cesses in p + A and A + A collisions. 

The initial energy density at RHIC is expected to 
be higher than at SPS. If one observes significant sup- 
pression of large-pT hadrons at RHIC as was predicted 
H [ui| , it clearly reveals an initial condition dramatically 
different from the CERN SPS. 

The observed high-p^ pion spectra in central Pb + Pb 
collisions cannot be due to collective hydrodynamic flow, 
since there will always be high-p^ partons produced in 
the coronal region of the two overlapped nuclei where jet 
propagation and fragmentation will not be influenced by 
the dense matter. To verify that these spectra are from 
jet production and fragmentation rather than from hy- 
drodynamic flow, one can measure the azimuthal particle 
correlation (selecting particles above a certain px) rela- 
tive to a triggered high-p^ particle as was proposed in 
Ref. |}0| . One should see a double-peak structure charac- 
teristic of a jet profile. One can use this method at even 
moderate pt (where there are still not many particles 
per event) to determine the contribution from semi-hard 
processes and the px range for which use of a thermal 
fire-ball model is justified. Otherwise, the extracted tem- 
perature and radial flow velocity can be misleading. 
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